Charge conservation in RHIC and contributiuons to local parity violation observables 
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Relativistic heavy ion collisions provide laboratory environments from which one can study the 
creation of a novel state of matter, the quark gluon plasma. The existence of such a state is postu- 
lated to alter the mechanism and evolution of charge production, which then becomes manifest in 
charge correlations. We study the separation of balancing charges at kinetic freeze-out by analyzing 
recent result on balancing charge correlations for Au+Au collisions at ^smn = 200 GeV. We find 
that balancing charges are emitted from significantly smaller regions in central collisions compared 
to peripheral collisions. The results indicate that charge diffusion is small and that the centrality 
dependence points to a change of the production mechanism. In addition we calculate the contri- 
butions from charge-balance correlations to STAR's local parity violation observable. We find that 
local charge conservation, when combined with elliptic flow, explains much of STAR's measurement. 



I. INTRODUCTION 



In relativistic heavy ion collisions a large fraction of the 
detected charge is produced throughout the time evolu- 
tion of the system. In particular for Au+Au collisions 
at y/s^N = 200 GeV the charged particle multiplicity is 
initially N c h = 158, whereas in the final state dN c h/drj 
is observed to be many hundreds of particles per unit 
rapidity, over multiple units of rapidity, with the precise 
values dependent on the impact parameter [1]. Despite 
the necessity of a correct treatment of the charge produc- 
tion inherent to relativistic heavy ion collisions, different 
reaction models employ very different production mech- 
anisms. In event generators like HIJING [2] and par- 
tonic/hadronic cascades such as URQMD [3] the charges 
are created early by the fragmentation of gluonic strings. 
In contrast, hydrodynamic models have the inherent as- 
sumption that a large fraction of the charge is produced 
at a later time at the deconfinement phase transition, i.e. 
when hadrons appear in the collision. Possible scenarios 
for the late stage production of charge include formation 
of hadrons from gluons, conversion of the nonperturba- 
tive vacuum energy into particles as well as quark pro- 
duction induced by the hadronization of quarks. In this 
context the late stage production of quarks has been pro- 
posed as a signature of the existence and lifetime of the 
quark gluon plasma 4 . In principle, charge production 
could occur at any time before the chemical freeze-out 
of the system. The problem is therefore related to the 
dynamics of the non-equilibrium phase, the quark-gluon 
plasma and finally the deconfinement phase transition. 
However, microscopic approaches from first principles are 
still in their infancy and can not yet provide further in- 
sight with regard to charge production mechanisms. It is 
therefore of great importance to explore the experimental 
possibilities to obtain valuable information. In the con- 
text of the search for late stage production of charges it 
has been proposed to investigate balancing charge corre- 



lations [4] . The idea is based on the principle that charge 
is created in balancing pairs that originate from the same 
point in space-time. In the subsequent rescattering of 
the charge carrier, which in principle can be hadronic 
or partonic, the balancing partners are then spread out 
within some finite distance of each other. However the 
motion is in general highly collective, i.e. the collision 
rates are high compared to the rate of expansion, so that 
the original correlation in space-time transforms into a 
correlation in momentum space in the emission profile 
at freeze-out. In a transport model the motion of the 
balancing partners could be separated into a collective 
mode, due the collective expansion of the system, and a 
diffusive relative mode, due to the collisions with other 
particles. The kinetic freeze-out configuration is then de- 
scribed by the single particle distributions with the ad- 
ditional constraint of having balancing partners emitted 
within a certain distance of each others, related to the dif- 
fusion constant and the available timescale for diffusion 
[3]. If charge is created late in the collision and diffusion 
is small the balancing partners are expected to be emit- 
ted within a small range and the correlation is strongest. 
The study of charge-balance correlations hence gives in- 
sight into the production and diffusion of charge. 

The observable that is employed to investigate these 
correlations is the charge balance function, which in the 
most general form is given by 
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Here dM/dpi denotes the differential charged particle 

multiplicity and AT_| (P2IP1) is the number of pairs with 

a positive charge emitted with momentum p\ and a neg- 
ative charge emitted with momentum pi and accordingly 
for N ++ , N |_ and N The balance function is con- 
structed in such a way that the like-sign subtraction sta- 
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tistically isolates the balancing partner. If we assume 
that all other correlations are charge insensitive, i.e. iden- 
tical for same-sign and like-sign pairs, then these corre- 
lations are removed from the balance function by sub- 
tracting the like-sign pairs from the opposite-sign pairs. 
In reality the additional sources of correlations are charge 
sensitive. In particular it has been shown that for small 
relative momenta the balance function is sensitive to cor- 
relations induced by the final-state interactions between 
the two particles due to Coulomb interactions and iden- 
tical particle symmetrization i5 . One needs to account 
for such correlations whenever p\ sa p 2 in Eq. above. 
The balance function in ([I]) is a six-dimensional function 
of the particles momenta. In the context of studies of the 
separation of balancing charges it is sensible to reduce the 
discussion to measures of the differences of the momenta, 
Pi~P2- In particular we will focus on the charge balance 
function in relative pseudorapidity Ar] and relative az- 
imuthal angle A<j>, which are given by 

B(A V ) = ~J dp 1 ~dp 2 B(p 1 \p 2 )6(Ar)-\r h - m \l?) 

and accordingly for A<j>. Balance functions of this form 
have been discussed in [?J [5j and been measured by the 
STAR collaboration [6l [7]. The experimental results 
have shown that balancing charges are in general highly 
correlated at freeze-out. The observed correlations are 
stronger for more central collisions compared to more 
peripheral collisions and remarkably it was found that 
for central collisions the observed charge-balance corre- 
lations are consistent with thermal emission of balancing 
pairs from the same collective velocity, i.e. perfectly lo- 
cal charge conservation at freeze-out [7] . However it 
could not be distinguished, how much of the observed 
narrowing of the balance function is due to a change in 
temperature and flow, as opposed to a smaller separation 
of balancing charges at freeze-out in coordinate space. 
In this paper we extend previous studies by performing a 
systematic analysis of the locality of charge conservation 
at break up based on the most recent STAR measure- 
ment of charge balance functions \J\. We modify STAR's 
single particle blast-wave model [8] to account for local 
charge conservation at freeze-out. In this model particles 
are emitted according to the blast-wave description with 
the additional constraint of local charge balance within 
a finite range of longitudinal rapidity a n and azimuthal 
angle a^. We then extract the separation of balancing 
charges by adjusting a v and to reproduce the observed 
charge balance functions in relative pseudo-rapidity Ar\ 
and relative azimuthal angle Acj>. Since the blast-wave 
parameters were adjusted by STAR to match spectra and 
elliptic flow for each centrality bin, this allows us to dis- 
tinguish between effects that are caused by a change in 
the freeze-out temperature and the collective flow versus 
effects that are due to the spatial separation of balancing 
charges, which is only affected by the production mech- 



anism and the subsequent charge diffusion. 
In addition, we show how charge-balance correlations 
contribute to STAR's local parity violation observable 
[S] and use the same model to perform a detailed anal- 
ysis of these "background" contributions. This paper is 
organized as follows: The modified blast-wave model is 
described in detail in Sec. II, results on charge separation 
are presented in Sec. III. In Sec. IV we give an intro- 
duction to the phenomenon of local parity violation in 
relativistic heavy ion collisions and show how the fluctu- 
ations of parity-odd observables measured at STAR are 
affected by charge-balance correlations. The findings are 
summarized in Sec. V. 

II. FREEZE-OUT MODEL 

In common parameterizations of the kinetic freeze-out 
configuration, the system is described by a local thermal 
equilibrium, 

where Xs( x ) is the characteristic function of the freeze- 
out surface S. The system is then completely described 
by the collective velocity profile u^ix), the kinetic freeze- 
out temperature T], in and the freeze-out surface S, which 
is usually chosen to be some volume at a fixed time in the 
lab frame. The conventional way to apply such a blast- 
wave model is then to choose a point on the freeze-out 
surface, create a particle with momentum according to 
the kinetic freeze-out temperature and the collective ve- 
locity, calculate the observable of interest and then repeat 
for several particles until sufficient statistics is achieved. 
The chemical composition is either governed by some pa- 
rameterization of the chemical freeze-out configuration or 
enforced to match experimental spectra. In many cases, 
e.g. when only single particle observables are consid- 
ered, this procedure is sufficient to phenomcnologically 
describe a variety of aspects of the hadronic spectra (see 
e.g. [5]). With regard to charge-balance correlations one 
needs to additionally incorporate local charge conserva- 
tion. This can be achieved in the following way: Instead 
of generating a single particle at a time, we generate an 
ensemble of particles with exactly conserved charges, for 
a given ensemble every particle is then assigned a col- 
lective velocity Vi such that all the follow the single 
particle blast-wave parameterization with the additional 
constraint of being emitted within a certain distance of 
each others. In the limit where charge conservation at 
kinetic freeze-out is perfectly local, this distance equals 
zero and all particles within a given ensemble are emit- 
ted as if their sources had the same collective velocity. 
If in contrast the balancing partner charges were spread 
out over the entire volume of the system, the locality 
constraint drops out and one recovers the conventional 



single particle blast-wave model. This extended freeze- 
out model then consists of three separate components, 
a parameterization of the chemical freeze-out configura- 
tion with exact charge conservation, a parameterization 
of the single particle freeze-out properties and a param- 
eterization of the separation of the particles within an 
ensemble. 

For the chemical composition we consider canonical en- 
sembles at a chemical freeze-out temperature of T c h em = 
175 MeV [TU] with total electric charge, strangeness and 
baryon number all equal to zero. We use a dilute descrip- 
tion of the hadron resonance gas to calculate the canon- 
ical partition functions according to the procedure out- 
lined in [5J. Finally we generate ensembles with zero net 
charge according to the canonical partition function. We 
consider mesons from the flavor octet and singlet ground 
state pseudoscalars and pseudovectors. Baryons are cho- 
sen from the ground state decuplet and octet. 
For the single particle freeze-out properties we use the 
STAR parameterization suggested in [H] that was shown 
to reproduce transverse momentum spectra as well as 
elliptic flow measurements. The transverse boost pro- 
file is given by an elliptic filled shell blast-wave, while 
the parameterization is boost invariant in the longitu- 
dinal direction. The parameters of the model are the 
kinetic freeze-out temperature 7fc in , the maximum trans- 
verse collective rapidities in-plane y x and out-of-plane y y 
and the in-plane and out-of-plane radii R x and R y of the 
freeze-out surface. The points of emission (x, y) are cho- 
sen uniformly within an ellipse with radii R x and R y . 
For a given point (x, y) on the freeze-out surface, the 
collective transverse rapidity is given by 



Vt = r (po + P2 cos{2(f> B )) , 




where we introduced po = (y x + 2/^)/2 and P2 — 
{Vx — Vz)l^ according to the notation in [8 . The 
direction of the transverse collective velocity is per- 
pendicular to the surface of the ellipse at the point 
(x,y), i.e. the collective velocity is given by it M = 
(cosh(y t ),sinh(y t )cos(0 B ),sinh(y t )sin((/)B),O). The sit- 
uation is illustrated in Fig. [I] We note that the parame- 
terization is only sensitive to the ratio of the in-plane to 
the out-of-plane radius R x /R y which has been extracted 
in [5], however the absolute values become meaningful 
when final-state interactions are considered. Hence this 
is addressed separately in Sec. III. a). 

To incorporate local charge conservation one has to 
generate collective velocities for all particles of an ensem- 
ble in such a way that they satisfy a specified relative 
distribution, while the properties of the single particle 
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FIG. 1. Visualization of the blast-wave parameterization in 
the transverse plane. At a point (x, y) in the transverse plane 
the associated collective velocity is perpendicular to the sur- 
face of the ellipse. According to Q the magnitude of the 
transverse collective rapidity increases linearly with the re- 
duced radius and has a second order harmonic component 
w.r.t the boost angle (f>B- The relative distribution of the bal- 
ancing charges |5| in the transverse plane is a Gaussian in 
A4>. 



distributions remain unchanged. For the relative distri- 
bution of the resonances within an ensemble we assume 
a gaussian distribution of the form 

-(-^rW-^f) < 5 » 

where (pi — atan \(R X /R y ) (yi/xi)] as shown in Fig. [l] 
Here Ay zi j denotes the difference in longitudinal rapid- 
ity between the particles i and j, and accordingly for 
A(f>ij. This particular choice has two important prop- 
erties that greatly simplify the sampling: a) The single 
particle distributions P(y z ,i) and P{4>i) are uniform in 
y z (longitudinal boost invariance) and (j) respectively, b) 
The use of a boost invariant blast-wave parameterization 
in combination with a relative distribution that factor- 
izes into a longitudinal and a transverse part allows us 
to generate the longitudinal and transverse collective ve- 
locities separately. 

For a given ensemble of N particles the collective veloc- 
ities can then efficiently be generated by a Metropolis 
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algorithm, with a target distribution 

[ da{[[P(y z ^ N e,J-^^) , (6) 
J i=l \ v / 

and accordingly for </>. It is straightforward to verify that 
this reproduces the correct single particle distributions 
in the limit of av,^ <C 1. In practice we verify empir- 
ically that this is the case for all calculations presented 
in the upcoming sections by checking transverse momen- 
tum and elliptic flow. The integration over a in Eq. © 
is part of the Metropolis routine, where the variable a 
effectively represents the average y z j and is introduced 
to accelerate convergence. 

Finally the resonances and particles are assigned their 
thermal momenta according to the kinetic freeze-out tem- 
perature and their respective collective velocity. Reso- 
nances are then decayed according to measured branch- 
ing ratios and lifetimes 11 . Since charge conservation 
is enforced on an ensemble-by-ensemble basis, particles 
originating from different ensembles are uncorrelated. 
Charge balance correlations can therefore efficiently be 
calculated on an ensemble-by-ensemble basis. 
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FIG. 2. Balance Function in A77 for different centralities. The 
STAR data is shown in blue lines, the red curve is obtained 
from the modified blast- wave model (see Sec. II). The balanc- 
ing charge separation in longitudinal rapidity a v is indicated 
for each centrality bin. The broadening of the balance func- 
tions for less central collisions is a result of both the higher 
kinetic freeze-out temperature and a larger spatial separation 
of balancing charges at freeze-out. 



III. SEPARATION OF BALANCING CHARGES 

The separation of balancing charges at freeze-out is ex- 
tracted from the experimental results [7J in the following 
way. We calculate balance functions for all charged parti- 
cles in relative pseudorapidity Arj and relative azimuthal 
angle A<p with the freeze-out model described in the pre- 
vious section. We then adjust the widths of the relative 
distributions a v and for each centrality bin, in order 
to match the STAR results. The single-particle blast- 
wave parameters are taken from STAR [5] , separately for 
each centrality. We assume perfect detector efficiency 
for the blast-wave calculation and use acceptance cuts of 
0.2 GcV/c < p t < 2.0 GeV/c and \r]\ < 1.0 as in the 
STAR measurement [7J. The necessary efficiency correc- 
tion is performed by normalizing all balance functions 
in such a way that they integrate to unity. In Fig. [2] 
we present the balance function in Arj compared to the 
STAR data for various centralities. Here 0% centrality 
corresponds to zero impact parameter and we refer to 
[7J for more details on the classification. The considered 
charge separations are indicated in the figure and corre- 
sponds to respective best x 2_nt - Fig. [3] shows the same 
results for balance functions in A<j). 

The charge-balance correlations in A77 show convinc- 
ing agreement with the STAR data. We find that the 
observed broadening of the balance function for more pe- 
ripheral collisions is a result of the higher kinetic freeze- 
out temperature and a larger separation of the balanc- 
ing charges. For more central collisions the STAR data 
show an additional peak at small relative pseudorapidi- 
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FIG. 3. Balance Function in Acj> for different centralities. 
The STAR data is shown in blue lines, the red curve is ob- 
tained from the modified blast-wave model (see Sec. II). The 
balancing charge separation in azimuthal angle cr^, is indi- 
cated for each centrality bin. The broadening of the balance 
functions for less central collisions is a result of the higher ki- 
netic freeze-out temperature, less transverse flow and a larger 
separation of balancing charges at freeze-out. 



ties which is not explained by the model. The overall be- 
havior is nevertheless very well described by only consid- 
ering charge-balance correlations. Also the consideration 
of additional sources of correlations has no significant 
effect on the result inferred for the longitudinal charge 
separation at freeze-out. 
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The agreement with the data is worse with regard to bal- 
ance functions in A<p. First of all the STAR data shows 
oscillatory behavior for the three most central bins. This 
is expected to be an artifact caused by sector boundaries 
of the detector as discussed in great detail in the experi- 
mental analysis [7|. However this issue complicates a re- 
liable extraction of the charge separation. We note that 
the behavior at large A</> is nevertheless well reproduced, 
even in this centrality region. Similar to the balance 
functions in Arj we find that the STAR data shows an 
additional peak at small relative angles for very central 
collisions again not predicted by the model. In addition 
there is a bump in the experimental data at small relative 
angles for very peripheral collisions. We will show in Sec. 
III. a) how this arises from the distortion of the balance 
function when correlations due to final-state interactions 
are taken into account. 

We emphasize that the change in the kinetic freeze-out 
temperature and the collective flow alone fail to explain 
the observed narrowing of the balance function for more 
central collisions. This feature of the STAR data can only 
be explained when the separation of balancing charges is 
taken into account. In Fig. [4] we present our result on 
the separation of balancing charges at freeze-out. The 
results are obtained from the fit to charge balance func- 
tions. The separation of balancing charges in longitudinal 
rapidity is smallest for central collisions and increases for 
higher impact parameters. The qualitative behavior for 
the separation of balancing charges in azimuthal angle 
is similar, however there is a bigger uncertainty in the 
results primarily arising from the ambiguities with the 
experimental data. We restrain from showing error bars, 
as they could only include statistical errors and discuss 
the systematic uncertainties associated with the calcula- 
tion instead. 

The longitudinal results arc fairly model independent, 
as the dependence on the parameterization in the trans- 
verse plane is weak and longitudinal boost invariance at 
mid-rapidity (recall that the analyzed data is obtained 
for |?7| < 1) is confirmed by various experimental obser- 
vations. In contrast we expect a significant model de- 
pendence of the results in azimuthal angle, as the em- 
ployed parameterization (see Eq. ([5])) is somewhat arbi- 
trary. Further insight with regard to this question might 
be achieved by the experimental study of more differ- 
ential observables such as two dimensional charge bal- 
ance functions in <fi and A(f>. Finally there is some un- 
certainty associated with the chemical composition and 
the proper treatment of decays, in particular related to 
DCA cuts and detector efficiency, which in principle can 
be addressed by incorporating a more detailed efficiency 
and acceptance model into the calculation. We expect 
those issues to affect our inferences of and a v on the 
order of 10-20% percent. 
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FIG. 4. Separation of balancing charges in longitudinal ra- 
pidity a v and azimuthal angle o<f, at freeze-out as a function 
vs. centrality of the collision. The balancing charges are least 
separated both longitudinally and transversely for central col- 
lisions. 



a) INFLUENCE OF ADDITIONAL SOURCES OF 
CORRELATIONS 

In the previous discussion of balance functions it has 
been assumed that the balancing partner charge is statis- 
tically isolated by the like-sign subtraction iV + _ — N ++ , 
i.e. only the balancing partner contributes to this expres- 
sion. This assumption does not hold when additional 
sources of correlations are taken into account and we 
have seen from the results presented in Fig. |3j that local 
charge conservation alone fails to explain the shape of the 
balance function for more peripheral collisions. Here we 
will consider final-state interactions between pion pairs 
to show that much of the experimentally observed dis- 
crepancy seen in Fig. [3] is a result of final-state interac- 
tions. In order to estimate the distortion of the balance 
function from residual interactions we follow the method 
provided in [5] . For every pair of balancing particles, one 
has to consider the correlation with other pairs of balanc- 
ing particles. The correlation weight for a 7r + 7r~ pair, p a 
and with another 7r + 7r~ pair, p c and pd, is then 

w(Pa,Pb;Pc,Pd) = \^ac{Pa,Pc)\ 2 \lpad(Pa,Pd)\ 2 

\Tpbc{PbiPc)\ 2 \^bd{PbiPd)\ 2 (7) 

where ip a p is the Coulomb wavefunction with asymptotic 
momentum p a — pp. For two pions of the same sign, the 
wavefunctions are symmetrized. The distortion of the 
balance function due to identical particle symmetrization 
and Coulomb correlations is then estimated from the ad- 
ditional contribution to the balance function numerator 
from the two pairs, where instead of increasing the pair 
distributions N ++ and in unit steps, the distribu- 
tions are increased by the respective correlation weight 
w(pa,Pb',Pc,Pd)- The pairs are generated according to 
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the blast-wave prescription described in Sec. II. We re- 
strict to 7r + 7r~ pairs and assume that the two particles 
from the same pair are emitted from the same collective 
velocity. The charged particle multiplicities are chosen 
according to STAR measurements [T] and we choose a 

source radius of yj R"i + Ry = 13 fm in the transverse 
plane. The longitudinal size of the system is determined 
by the freeze-out time which is chosen to be r = 10 fm/c 
to reproduce HBT spectra [5] ■ The distortion of the bal- 
ance is then estimated by applying the corrections to 25% 
of all pairs, in order to account for pions emitted from 
long-lived decays and because only a fraction of pairs in- 
volve two pions. This estimate is rather crude, but should 
be sufficient to see if the effects of final-state interactions 
are of the appropriate magnitude and shape to describe 
the discrepancies seen in Fig. [3j 

In Fig. [5] we show how the balance function for 50- 
60% centrality is distorted due to symmetrization and 
Coulomb correlations of ir + ir~ pairs. We find that the dip 
at small A0 present in the STAR data, can be reproduced 
when identical particle symmetrization and Coulomb cor- 
relations are taken into account. However this happens 
at the expense of higher values of the balance function 
at large A(f>, and therefore also affects the separation of 
balancing charges extracted in the previous section. In 
particular for the 50-60% centrality bin we find that the 
distortion of the balance function lowers ct^ by about ten 
percent. The effects of additional sources of correlations 
on the balance function in Ai] are expected to be signifi- 
cantly lower. This is because pairs of particles with small 
A<j) are also within two units of pscudorapidity, whereas 
for a pair of particles with small Arj they the particles 
can be spread out over the entire angle of 2tt. Hence one 
would naturally expect the effects to be smaller by a fac- 
tor of 7r. For a more detailed study of the residual sources 
of correlations non-resonant contributions of the strong 
interaction also have to be taken into account and the 
analysis has to be extended to other particle species. We 
recommend performing such studies in terms of balance 
functions of the invariant momentum of the pair q.i nv and 
its components q out , q slde and 



IV. CONTRIBUTIONS TO LOCAL PARITY 
VIOLATION OBSERVABLES 

The prospect of observing a signature of topologically 
non-trivial QCD gauge field configurations in relativistic 
heavy ion collisions has recently gained great attention 
H2] . It was shown that in the presence of such config- 
urations P- and CP- violating processes occur locally in 
regions with non-zero QCD topological charge [15] . which 
are expected to be created in the early stage of relativis- 
tic heavy ion collisions. In particular in the presence of a 
magnetic field positive charge is separated from negative 
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FIG. 5. Balance function B(A<f>) for 50-60% centrality from 
the blast- wave model (red curve). The corrections due to 
final-state interactions are included in the blue curve, STAR 
data is shown in black. The distortion of the balance function 
due to symmetrization and Coulomb interactions explains the 
behavior of the STAR data at small relative angles. 

charge along the magnetic field [14] . This phenomenon is 
referred to as the "Chiral magnetic effect" . In the situa- 
tion of relativistic heavy ion collisions the magnetic field 
is induced by the colliding nuclei. This implies that par- 
ity is locally violated on an event by event basis as pos- 
itive and negative charges separate along the magnetic 
field, which is perpendicular to the reaction plane (out- 
of-plane) [S]. However because there is no direct V- and 
CP- violation in QCD, the probability to generate charge 
separation parallel or anti-parallel to the magnetic field 
is equal and the expectation value of any P-odd observ- 
able vanishes. The possibilities to identify the existence 
of a local parity violating effect are therefore confined to 
the study of fluctuations of "P-odd observables (which are 
then P-even). The observable measured by STAR is 

7a '" = ^ Mjtp ' (8) 

where a and /3 represent positive or negative charge, 
M a and Mp are the corresponding multiplicities, the 
azimuthal angles <p are measured about the beam axis 
and ipRp is the angle of the reaction plane which we 
will set to zero in the following without loss of gener- 
ality. The observable is "P-even, so that in principle 
other sources of correlations can contribute to the sig- 
nal. Since one is looking for charge separation along 
the out-of-plane direction, the idea is to compare the 
charged particle correlations out-of-plane to the corre- 
sponding in-plane-correlations. By use of the identity 
cos(4>i + (j)j) = cos(0i) cos(0j) — sin(0j) sin(<pj), the STAR 
observable subtracts in-plane from out-of-plane correla- 
tions to eliminate other sources of correlation that are in- 
dependent of the orientation w.r.t the reaction plane. If 
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there are no reaction-plane-dependent background con- 
tributions, charge separation in out-of-plane direction 
then causes 7 to be negative for same-sign pairs and pos- 
itive for opposite-sign pairs. The correlators in Eq. ^ 
have recently been measured by the STAR collaboration 
[|5]. In addition to the integrated signal differential mea- 
surements have been performed in relative pseudorapid- 
ity Arj and transverse momentum p t |9j. The reported 
results are qualitatively in line with various expectations 
from the chiral magnetic effect and significantly larger 
than reaction-plane-dependent background contributions 
from MEVSIM, HIJING and UrQMD simulations 0. 
However there is an ongoing discussion on the interpre- 
tation of the data, in particular if the chiral magnetic 
effect can produce a signal of the observed order of mag- 
nitude [T5J [TS] and whether there are more "traditional" 
explanations of the observed signal (e.g. cluster parti- 
cle correlations [17]; charge conservation and flow [18], 
momentum conservation [16 [ I19j). Furthermore the sep- 
arate analysis of in-plane and out-of-plane correlations 
has revealed, that the correlations for same-sign particles 
are mainly in-plane and back-to-back I20j , indicating the 
presence of other sources of correlations, and additional 
measurements have been proposed to clarify the situation 
[2T] , In the following we will show how charge balance 
gives rise to reaction-plane-dependent correlations and 
calculate the respective contributions to the local parity 
violation observable from the blast- wave model described 
in Sec. II. 

From the previous discussion (see Sec. I) it is clear that 
charge-balance correlations only contribute to the differ- 
ence of opposite-sign and same-sign pair correlations. We 
will therefore only discuss contributions to 
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IP 



(27+- -7++ -7—)' 



(9) 



The observable 7p then compares in-plane vs. out-plane 
correlations for opposite-sign vs. same-sign charged 
pairs. In this context the discussion of charge-balance 
correlations has to be extended to balance functions sen- 
sitive to the angle of the pair with respect to the reaction 
plane cf>, i.e. the object of interest is 



B(0,A0) 



1 



dM/d<t> 



dv\^r—dp2 B(pi,p 2 ) (10) 
dpi 



Balance functions B(<fi, A<p) are presented as a function of 
A(j) in Fig. [6] for events with centralities of 40-50%. The 
results are obtained from the blast-wave model described 
in Sec. II for the charge-separation parameters extracted 
in Sec. III. The balance function for <j> = 0° (in-plane) 
is narrower than the balance function for = 90° (out- 
of-plane). The stronger focusing of balancing charges 
derives from the greater collective flow in-plane vs. out- 
of-plane. For = 45°, the distribution is biased toward 
negative values of A(j>. This is expected given the elliptic 



asymmetry, vi > 0, which leads to more balancing parti- 
cles toward the — 0° direction as opposed to = 90°. 
Depending on which quadrant is located, the balanc- 
ing charge tends to be found more toward = 0° or 
(6 = 180°. 



0.25 




180° 



FIG. 6. Balance function B(0,A0) for 40-50% centrality 
shown as function of the relative angle included by balancing 
partners for = 0° (black squares), 45° (red triangles) and 
90° (blue circles). The balance function is narrower for in- 
plane pairs compared to out-of-plane pairs. For intermediate 
angles the balance function is biased toward negative angles. 

The width and asymmetry of the balance function can 
be quantified by the moments, 

c 6 (0) = !tt f dA<j> 5(0, A0) cos(A0), (f f ) 
ZbW J 



where 



s b {<i>) = / dA0 5(0, A0)sin(A0), 



z b {<t>)= / dA0 5(0,A0), 



(12) 



is the normalization of the balance function and repre- 
sents the probability of detecting the balancing charge 
given the observation of a charge at 0. It would be 
unity for a perfect detector, but is reduced by both 
the finite acceptance and efficiency of the experiment. 
The moments Cb(0) and s;,(0) are displayed in Fig. 
[7] The quantity Cf,(0) expresses the width of the 
balance function and would be unity for a very narrow 
balance function whereas it vanishes in the case where 
the balancing charges were emitted randomly. The 
quantity Sb(0) measures the degree to which the balance 
function is asymmetric under a reflection symmetry of 
A0 — > — Acf>. For pairs around = 45° this corresponds 
to the probability for the balancing charge to be emitted 
in in-plane direction vs. in out-of-plane direction. 

In order to relate the correlation jp to moments of the 
balance function 5(0, A0), some remarks on the event 
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FIG. 7. The moments of the balance function, Cb(c/>) and 
Sb{(j>), represent averages of cos(A^)) and sin(A0) across the 
balance function. These are plotted as a function of <f> for vari- 
ous centralities. The structure of Cb(<fi), which is maximized at 
cj> — 0° , 180°, illustrates how the balance function is narrower 
for in-plane emission and for more central collisions, while the 
structure of Sb(4»), which is positive for (f> = 135°, 315° and 
negative for (j> — 45°, 225°, shows how balancing charges pre- 
fer to be emitted in the in-plane direction. The oscillations 
increase for more peripheral collisions. 



averaging are in order. While the event average of 7p 
is simply given by (jp) [9 |; charge balance functions are 
calculated as the ratio of event averages, 



B(<p,A<p) 



(jV + _(0|A0)-A ++ (0|A0)) 
(dM/d(f>) 
(iV- + MA0)-iV— QflAfl) 
{dM/dcj)) 



(13) 



By use of the angle addition formula cos(0,; + <f>j) = 
cos(20i) cos(A0) — sm(2cf>i) sin(A</>) we find the expres- 
sion relating the correlation "fp to moments of the bal- 
ance function, 



<MV> 
(M) 



2 

W) 



dA<t) 



cos(2</>) cos(A0) 



(14) 

sin(2(^>) sin(A<^)] , 



where it has been assumed that there are equal numbers 
of positive and negative charges and the factor of 1/ (M) 
makes the result independent of the multiplicity. In the 
limit of very small multiplicity bins, the left-hand side of 
Eq. 14 could be replaced by (M) (jp) 
It is insightful to express the correlations, jp , in terms 
of the moments of the balance function defined in Eq. 



(HI, 



(M 2 7P ) 
2(M) 



Vl{Cb{4>)) + v 2c - v 2s 



(15) 



v 2c 

V2s 



(/(</>)> = 



{c b {4>) cos(20)) 
(s b {<t>) sin(20)), 
1 f , , dM 



v 2 {c b {(j>)). 



(M) 



(16) 



The three contributions to 7p derive from: a) having 
more balancing pairs in-plane than out-of-plane (v 2 (c b )), 
b) having the in-plane pairs being more tightly correlated 
in A(j) than the out-of-plane pairs (v 2c ) and c) having the 
balancing charge more likely being emitted toward the 
event plane (v 2s ). 

The three contributions to the signal obtained from the 
blast-wave calculation are displayed in Fig. [8| We as- 
sume perfect detector efficiency for the blast-wave calcu- 
lation and use the same acceptance cuts in transverse mo- 
mentum and pseudorapidity, i.e. 0.15 GeV < p t < 2 GeV 
and I r\ \ < 1. The necessary efficiency correction is done 
by rescaling the results to reproduce the experimental 
normalization of the balance function [7], i.e. we multi- 
ply the expressions for (c b ((f))), v 2c and v 2s by the ratio 
of experimental to blast-wave normalization. In addition 
to the results obtained for the (realistic) charge separa- 
tion extracted in Sec. Ill, we present the contributions 
assuming that charge conservation is perfectly local at 
freeze-out. These are the strongest possible contributions 
to the signal and should give an upper limit. In order to 
compare to STAR data [5] we modify the left hand side 
of ( 15 1 in the following way. Since Af 7p is independent of 
the multiplicity we expect (M 2 7p) w (M)(M^p), with 
< M7P > contained in M (jp) within errorbars. The left 
hand side of ( 15 1 then simplifies to M/2('jp), where M is 
the experimental multiplicity for a single event account- 
ing for efficiency and acceptance of the detector [2"2"| . 

The contributions of charge-balance correlations for a 
realistic charge separation are of the same size as the ex- 
perimental signal and exhibit similar qualitative behavior 
with respect to the centrality dependence. The system- 
atic error associated with this prediction originates pre- 
dominantly from the particular parameterization of the 
charge separation in azimuthal angle (see Sec. II, Eq. 
^ for more details) and we expect the conclusions to be 
reliable on a ten percent accuracy level. 
In the present discussion of the interpretation of STAR 
results it has been proposed to analyze the in-plane and 
out-of-plane correlations separately [20 . With respect 
to the difference between opposite-sign and same-sign 
correlations this becomes trivial as the sum of in-plane 
and out-of-plane correlations is related to reaction-plane- 
independent balance functions by 



- <t>j))+- - (C0s((f>i - (j)j))- 

dAcj) B(A<j>) cos(A0) , 



(17) 



which are well understood from the physics of charge 
balance as discussed in detail in the previous section. 
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FIG. 8. Difference between opposite-sign and same-sign par- 
ity observable from STAR (black dots) and blast-wave calcu- 
lations for realistic charge separation at freeze-out (red dots) 
and perfectly local charge conservation (blue dots). The three 
contributions to the signal are defined in Eq. ( 15 I and are 



plotted with dashed lines. V2{ct) (squares) derives from hav- 
ing more balancing pairs in-plane than out-of-plane while V2,c 
(triangles up) quantifies the degree to which in-plane pairs 
are more tightly correlated than out-of-plane pairs. V2, s (tri- 
angles down) reflects that the balancing charge is more likely 
to be found toward the event plane. 



The correlations can then be separated into their in-plane 
and out-of-plane components by the use of trigonometric 
identities |20j . both showing similar agreement with the 
STAR data [5] as the results presented in Fig. [8] 
In addition to the integrated signal the correlations have 
been analyzed differentially in pseudorapidity Arj and 
transverse momentum Pt,i+Pt,j [H]- In order to relate the 
differential correlators to moments of the balance func- 
tion, one has to account for the number of pairs in each 
differential bin. Hence Eq. 



(N P (A V ) 1P (A V )) 
(M) 



(M) 



D 



(|14| has to be modified to, 

dA<jy (18) 
, A<f>, Arj) cos(20 + A0) , 



-(f) 



where N p (Ar])/M 2 is the fraction of charged particle 
pairs in the respective pseudorapidity bin. The right 
hand side of |l9| ) can be obtained straightforward from 
the blast-wave model. In order to compare to experimen- 
tal data we will assume that, for a given event, the cor- 
relation 7p(A?7) scales with the total multiplicity rather 
than the fraction of pairs in the respective bin, i.e. 



M 2 



M 2 ^p(Ari) 



N P (Ar,) 

M 2 
NpjArj) 

M 2 



(MV(At7)) 



M(M){ 1P {A V )) 



where the average fraction of pairs per pseudorapidity 
bin can be obtained from the single particle spectra of 
the blast-wave model. The charge balance contributions 



to the differential observables are presented in Fig. [9] and 
10 respectively. The blast-wave results are obtained for 
realistic charge separation at freeze-out, the qualitative 
behavior is similar to what is observed at STAR. When 
regarded differentially in pseudorapidity the correlations 
are confined to a finite range in Arj. This is because 
balancing charges are likely to be emitted within a narrow 
range of rapidity. This can already be seen from balance 
functions in relative pseudorapidity (see Sec. II). When 
regarded differentially in p t the correlations increase with 
the momentum of the pair due to larger anisotropy (i.e. 
higher v?) and more collective flow. 
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FIG. 9. Differential Parity Observable from STAR for 30- 
50% centrality (black solid) and blast-wave calculations for 
30-40% (red dashed) and 40-50% centrality (blue dashed). 
The correlations decay in Arj as balancing charges tend to be 
emitted in a narrow range of relative pseudorapidity. 



V. CONCLUSIONS 

We have seen in Sec. Ill that the charge balance func- 
tions recently observed by the STAR collaboration [7] can 
be understood as a signal of highly localized charge bal- 
ance at kinetic freeze-out. The observed narrowing of the 
charge balance function for more central collisions, that 
could not be observed in URQMD and HIJING simu- 
lations 7J, is explained by a smaller separation of bal- 
ancing charges at freeze-out, which is not inherent to 
these models. The notion, that the observed narrow- 
ing of the balance function for more central collisions is 
entirely due to a reduced freeze-out temperature and in- 
creased transverse flow, can be dismissed. Even though a 
reduced freeze-out temperature and increased transverse 
flow lead to higher correlations of balancing charges in 
more central collisions, the relevant parameters describ- 
ing the flow and temperature of the breakup stage are 
constrained by spectra and elliptic flow measurements. 
As seen in Fig. [4j the separation of balancing charges in 
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FIG. 10. Differential Parity Observable from STAR for 30- 
50% centrality (black solid) and blast-wave calculations for 
30-40% (red dashed) and 40-50% centrality (blue dashed). 
The correlations increase in the transverse momentum of the 
pair p t = (pt,i +pt,j)/2, because of higher anisotropy and 
more collective flow. 

A77 appears smaller for central collisions than for periph- 
eral collisions. This separation is related to the separa- 
tion in coordinate space and the production time through 
the relation er J( 1=3 <t 2 /t. In p+p collisions the separation 
is driven by the dynamics of breaking gluon-strings, or 
color flux tubes. By tunneling out of the vacuum, balanc- 
ing charges typically appear separated by ~ 1/2 fm along 
the z-direction. If the particles are produced ~ 1/2 fm/c 
after the initial impact, they are already separated by a 
good fraction of a unit in pseudorapidity. This separa- 
tion can only increase in the subsequent time evolution, 
as the charges diffuse relative to one another. 
The observation of smaller a v for central Au+Au colli- 
sions implies that either: 

a) Charges are produced at later times, as this would 
allow the same a z to give a smaller a^. Late charge 
production is expected from delayed hadroniza- 
tion, as the formation processes of hadrons from 
quarks and gluons involve the creation of additional 
charge. 

b) Charges are produced early by a different mech- 
anism than the decay of color flux tube. In 
this situation cr z would be smaller initially, while 
r remains small. For instance, if the matter 
became isotropized gluonic matter at very early 
times, quark production might result from colli- 
sions rather than from tunneling through longitu- 
dinal flux tubes. This scenario also requires charge 
diffusion to be small, to constrain a n at freeze-out 
close to the initial separation at creation. 

Both of these explanations, or a combination of the two, 
might explain the smaller separation of balancing charges 



a n for more central collisions. The results presented 
in Sec. Ill therefore point to a change in the charge- 
production mechanism from p+p to central Au+Au col- 
lisions. However to draw further conclusions it is nec- 
essary to compare the extracted separations of balanc- 
ing charges to transport properties of the hadron reso- 
nance gas and the quark gluon plasma. As a first step 
it would be important to know how much the particles 
diffuse between chemical and kinetic freeze-out. These 
results could in principle be obtained from hadronic cas- 
cade models or transport theory. 

The separation of balancing charges in azimuthal angle, 
described by a^, shows a similar behavior with respect 
to the centrality dependence. As the size of the sys- 
tem in the transverse plane increases with centrality, this 
is qualitatively expected as long as the charges can not 
diffuse far away from one another, which would require 
both early production time and large diffusion. Hence 
the analysis of charge separation in azimuthal angle does 
not give any further insight concerning charge produc- 
tion mechanisms. 

With respect to STAR's local parity violation observ- 
able, we have seen that local charge conservation gives 
rise to a reaction-plane-dependent "background" signal, 
contributing to the difference of opposite-sign and same- 
sign correlations. The balancing charge correlations are 
of the same size as the observed correlations and exhibit 
similar qualitative behavior with respect to all measured 
properties. This suggests that local charge conservation 
is the dominant effect with regard to the difference of 
same-sign and opposite-sign parity observable. We em- 
phasize that this does not explain the experimentally ob- 
served strong same-sign correlations. However, it is clear 
that neither charge balance or parity fluctuations can 
produce same-sign correlations without affecting the dif- 
ference of opposite-sign and same-sign correlations. A 
mechanism that has been suggested to produce strong 
same-sign correlations is momentum conservation |16j . 
This is the topic of a separate study. 

This work was supported by the U.S. Department of 
Energy, Grant No. DE-FG02-03ER41259. 



[1] B.I. Abeleve et al. [STAR Collaboration], Phys. Rev. C 

79, 34909 (2009) 
[2] X.N. Wang and M. Gyulassy [HIJING], Phys. Rev. D 44, 

3501 (1991) 

[3] S.A.Bass et al [URQMD], Progr. Part. Nucl. Physics Vol. 
41, 225 (1998) 

[4] S.A. Bass, P. Danielewicz and S. Pratt, Phys. Rev. Lett. 

85, 2689 (2000). 
[5] S. Cheng, S. Petriconi, S. Pratt, M. Skoby, C. Gale, S. 

Jeon, V. Topor Pop and Q.-H. Zhang, Phys. Rev. C69, 

054906 (2004). 

[6] J. Adams et al. [STAR collaboration], Phys. Rev. Lett. 
90, 172301 (2003). 



11 



[7] M.M. Agg arwal et al. [STAR Collaboration], 

|arXiv:100 5.2307 [nucl-ex] (2010). 
[8] J. Adams et al. [STAR collaboration], Phys. Rev. C 72, 

14904 (2005). 

[9] B.I. Abelev et al. [STAR collaboration], Phys. Rev. Lett. 
103, 251601 (2009). 
[10] P. Braun-Munzinger et al., Phys. Lett. B 518, 41 (2001). 
[11] K. Hagiwara et al. [Particle Data Group], Phys. Rev. D 

66, 010001 (2002) 
[12] S.A. Voloshin, Phys. Rev. C 70, 057901 (2004) 
[13] D.E. Kharzeev, Nucl. Phys. A 830, 534C (2009) 



[14] D.E. Kharzeev, Annals Phys. 325, 205 (2010) 

[15] M. Asakawa, A. Majumder and B. Muller, 

|arXiv:100 3.2436 [hep-ph] (2010) 
[16] S. Pratt, |arXiv: 1002 1758V 1 [nucl-th] (2010). 
[17] F. Wang, TarXTv:0911. 1482^2 [nucl-ex] (2010) 
[18] S. Schlichting and S. Pratt, |arXTv:1005.5341| [nu cl-th] 
[19] A. Bzdak, V. Koch and J. Liao, arXiv: 1008.4919 [nucl-th] 
[20] A. Bzdak, V. Koch and J. Liao, Phys. Rev. C 81 031901, 

(2010) 

[21] S.A. Voloshin, [arXiv:1006.1020V l [nucl-th] (2010) 
[22] Experimental multiplicities were generously provided by 
the STAR collaboration. 



